Four-spin-exchange- and magnetic-field-induced chiral order in two-leg spin ladders 
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We propose a mechanism of a vector chiral long-range order in two-leg spin-^ and spin-1 antiferro- 
magnetic ladders with four-spin exchanges and a Zeeman term. It is known that for one-dimensional 
quantum systems, spontaneous breakdown of continuous symmetries is generally forbidden. Any 
vector chiral order hence does not appear in spin-rotationally [SU(2)]-symmetric spin ladders. How- 
ever, if a magnetic field is added along the S z axis of ladders and the SU(2) symmetry is reduced 
to the U(l) one, the z component of a vector chiral order can emerge with the remaining U(l) sym- 
metry unbroken. Making use of Abelian bosonization techniques, we actually show that a certain 
type of four-spin exchange can yield a vector chiral long-range order in spin-i and spin-1 ladders 
under a magnetic field. In the chiral-ordered phase, the Z2 interchain-parity (i.e., chain-exchange) 
symmetry is spontaneously broken. We also consider effects of perturbations breaking the parity 
symmetry. 

PACS numbers: 75.10. Jm,75.30.Kz,75.40.Cx,75. 50. Ee 



I. INTRODUCTION 



Quantum spin ladder systems have attracted much at- 
tention for more than a decade. 1 There are numerous 
compounds with a ladder structure. It is well known now 
that ladders contain considerably rich physics: for exam- 
ple, they provide fascinating ordered states, 2-5 string or- 
der parameters detecting different topological sectors, 6 
interesting exactly solvable points, 7 ' 8 etc. In particu- 
lar, physics of two-leg ladders has been greatly devel- 
oped. Low-energy properties of nonfrustrated two-leg 
ladders have been well understood: as soon as an in- 
finitesimal antiferromagnetic (AF) [ferromagnetic (FM)] 
interchain exchange is added between two spin-i AF crit- 
ical chains, a gap opens and a rung-singlet [Haldane] 
type spin-liquid state occurs, 6,9 while even if an inter- 
chain coupling with an arbitrary strength is introduced 
between two spin-1 AF Haldane-gapped chains, the gap 
always survives. 10 Based on these properties of standard 
ladders, many physicists recently have been focusing on 
more realistic or frustrated ladder systems. Among them, 
effects of four-spin exchanges have been studied inten- 
sively now. Although it is generally difficult to under- 
stand effects of four-spin exchanges compared with con- 
ventional two-spin ones, thanks to the discovery of a du- 
ality relation in two-leg spin-i ladders, 3 powerful ana- 
lytical tools (conformal field theory, Bethe ansatz, etc.) 
and accurate numerical methods in (1+1) dimensions, 
physical properties of four-spin exchanges in ladders have 
been gradually elucidated. 11 For instance, it has been 
shown that in spin-rotationally [SU(2)]-symmetric spin- 
5 ladders, four-spin exchanges can cause novel phases: 
a staggered dimer order, a scalar chiral one, and a spin 
liquid with a strong (but short-range) vector chirality 
correlation. 2-4 Moreover, in SU(2)-symmctric spin-1 lad- 
ders with a biquadratic term, a certain parameter re- 
gion where a spin-nematic correlation is dominant is pre- 



dicted. 12 

In this paper, motivated by the above studies of spin 
ladders, we investigate two-leg spin-^ and spin-1 AF lad- 
ders with four-spin exchanges in the presence of a mag- 
netic field, which reduces the SU(2) symmetry to a U(l) 
one. For these ladders, we propose a mechanism of a vec- 
tor chiral long-range order which spontaneously breaks 
the parity symmetry for the interchain (rung) direction. 

The paper is organized as follows. Throughout the pa- 
per, we apply the effective field-theory approach and fo- 
cus on the weak two- and four-spin rung-coupling regime. 
In Sec. II, we first review simple spin-i ladders without 
four-spin exchanges under a magnetic field using Abelian 
bosonization techniques and discuss the physical mean- 
ing of locking the bosonic fields. Then, based on the 
discussion, we prove the mechanism of a vector chiral 
order in the next two sections: Sees. Ill and IV are de- 
voted to investigating spin-i anc i spin-1 cases, respec- 
tively. In Sec. V, we shortly consider a few perturbation 
terms breaking the rung-parity symmetry: celebrated 
Dzyaloshinsky-Moriya (DM) interactions 13 and a parity- 
breaking three-spin term. 14,15 Finally, we summarize our 
results and then discuss them in Sec. VI. 



II. BOSONIZATION FOR SPIN-i LADDERS 



1,11 



We begin with a fundamental spin-i AF ladder, 
whose Hamiltonian is 



- H /]Sn.i< (!) 



where S n j is the spin-i operator on site j of the nth 
chain (n = 1,2), J > (Jj_) is the intrachain (rung) 
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coupling, and H is a magnetic field. At the decoupled 
point of J± = 0, the low-energy effective Hamiltonian for 
the nth chain is a Gaussian model 

nf = j dx v -[K-\d x <t> n f + K(d x e n )% 

where x = ja>o (ciq is the lattice constant), 4>n is the scalar 
boson field, and 8 n is the dual to 4> n : these two fields sat- 
isfy a canonical commutation relation [<j> n (x), d y 9 n {y)] = 
iS(x-y). Quantities v and K denote the spin- wave veloc- 
ity and the Tomonaga-Luttinger liquid (TLL) parameter, 
respectively. The value of K (v) runs from 1/2 (irJao/2) 
to 1 (0) with H increasing from to the upper critical 
field 2 J. Spin operators are also bosonized as 



S n,3 



M + ^=d x <t> n 

+(-1 YA 1 sin(\/4^„ + 27rMj) + • • • , (2a) 
+B 1 sin(V4^0„ + 2nMj) + •••], (2b) 



where M(H) — (S*j), and A q and B q are nonuniversal 
constants. 16 From this bosonization framework, the low- 
energy theory of the ladder (1) with | J±\ <C J is written 
as 

Weff = / [Kt\d x <t> t f + K t (d x 8t) 2 ] 
J t=± 1 



J± 



r A 2 

\Bl cos(V270_) + -± cos(v / 8tt</>-) 



cos(V87T0 + + AnMj) 



+ 



+ ■ 



MJ ± d x <p + 
(3) 



where we introduced new fields <j>± — (<p\ ± 4>i)j\[2 and 
6± = (6>i±6» 2 )/\/2 that obey [<j>±{x),d v e±{y)\ =iS(x-y). 
New TLL parameters K± and velocities v± are estimated 
as 



K± w k(\^K 



v± 



< 



1±K 



2ttv 
J±ap \ 
2wv )' 



(4) 



The boson linear term d x <p + can be absorbed into the 
Gaussian (boson quadratic) part by the shift <j) + — 

<j) + + a/ -K + M which accompanies the correction of 



the magnetization, (S n>j ) = {l-^^-K + )M = M. Here- 
after, we will again use symbols 4> + and M for 4> + and 
M, respectively. 

The U(l) spin rotation around the S z axis, the trans- 
lation by one site for the chain direction, the site-parity 
operation for the same direction and the rung-parity 
transformation (i.e., exchange between two chains) cor- 
respond to 6+ -> 6 + + const, {<j> + , 0+) -> (4> + + V2ir(M + 



i/2),e+ + v^F), (4>+(x),4>-{x),0±(x)) -> (-4>+(-x) + 

yftj2,-4>-(-x),6 ± (-x)) and (<£_,0_) - -(0_,0_)> 
respectively. The Hamiltonian (3) hence docs not con- 
tain 6 + vertex terms, <fi+ ones without spatially oscillat- 
ing factors e ± m47rM i ; anc j s j ne terms with <f>± and 0_. 

The form of Eq. (3) tells us that the low-energy physics 
of the ladder (1) can be described by two parts, (</>+, 0+) 
and (<fi_,6-) sectors. In the present notation, scaling 



are, respectively, 



dimensions of e m ^ 9t and e lm ^^ 
n 2 /(2K t ) and 2m 2 K t . The (<f)-,0-) sector always takes 
a massive spectrum due to relevant cosine terms, irre- 
spective of the value of M. Since cos(v / 27r#_) is usually 
more relevant than cos(\/87r0-) in the ladder, namely 
K- > 1/2, 17 the dual field 6- is pinned. For the AF- 
rung case, = ±y/ir/2, whereas for the FM-rung case, 
(#_) = 0. On the other hand, physics of the (</>+, 0+) sec- 
tor depends on M. When M — 0, i.e., H is smaller than 
the lower critical field Hi of the ladder, 18 cos(v / 87r0+) 
is relevant and a spin gap thus emerges. This state is 
nothing but the Haldane or rung-singlet spin liquids. It 
is known that pinning the field <f> + corresponds to the 
emergence of nonlocal string orders. 6 Inversely, if M ^ 
(H > Hi), cos(\/87r0+ + 47rMj) is irrelevant due to the 
factor e ±l4nMj and the gap vanishes. Therefore in the 
case of M ^ 0, the ((j> + ,6+) sector is described by a 
Gaussian model, which indeed corresponds to the field- 
induced TLL phase. These scenarios in the ladder (I) 
would be robust against small perturbations conserving 
symmetries of the ladder (e.g., XX Z anisotropy, further- 
neighbor exchanges, etc.). 

As typical order parameters for two-leg ladder systems, 
one can find three quantities: the z component of spin 
vector chirality Vj — (S\j x S2j) z and two magnetic 
moments Af±j — Sf j ± S^j- It is important to note that 
Vj and AT- j- are odd for the rung-parity operation, while 
Sf+,i is even. In the SU(2)-symmetric case with H = 0, 
which corresponds to the gapped spin-liquid state, ex- 
pectation values of these order parameters must be zero 
because the presence of them means the spontaneous 
breakdown of the continuous SU(2) symmetry and such 
a symmetry breaking is generally forbidden by Mermin- 
Wagner-Hohenberg argument. 19 The gapped spin- liquid 
state still remains if H < Hi and M = 0. Meanwhile, 
in the U(l)-symmetric case with M/0 their finite ex- 
pectation values can be allowed since they do not violate 
the U(l) symmetry. Formula (2) enables us to represent 
these order parameters as 



Vj w --Bjfsin(V27r0_) 



(5a) 



Af- 



(-iy2A 1 cos(V27T (? !) + + 2-kMj) sin^Tr^ 



+a \j-d x 4>- + 



(5b) 



JV+j ~ (-iy ? 2^isin(V27n/>+ + 2ttMj) cos(V2tt</> 



+2M + a \ -d x <j) + + 



(5c) 



We see that Vj and Af-j are actually odd for the rung- 
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parity transformation (</>_,#_) — > — ((/>_,#_), whereas 
A/+ j is even. The leading part of A/± consists of two 
fields <f>±. As we mentioned above, since </>+ is not locked 
in the case of M ^ 0, a magnetic order with (A/±) 7^ 
is shown to disappear in the case. On the other hand, 
remarkably, the leading term of Vj contains only the field 
that can be locked even in the U(l)-symmetric case 
with M ^ 0. 

From the above discussion on order parameters, it is 
inferred that the vector chiral order emerges in a certain 
class of U(l)-symmetric ladders under a magnetic field. 
For the case of the simple ladder (1), we find (Vj) = 
because of = or ±yJn/2. However, in other words, 
it is expected that if the 0_ -locked position is moved 
even slightly by additional perturbations, the chiral order 
would appear. We show below that some types of four- 
spin exchanges are indeed such a desirable perturbation. 

Here, we again remark on Eq. (5). As one saw above, 
for the simple spin-i ladder (1) with M = 0, boson fields 
in both (<j) + ,6+) and (0_,0_) sectors are locked. The 
physical meanings of locking </>_ or 9-, however, have not 
been discussed well so far. Equation (5) clearly shows 
that locking 9- can induce the vector chirality, while 
locking </>_ can do magnetic orders. Moreover, one finds 
that Eq. (5a) is analogous to the celebrated supercurrent 
formula in a Josephson junction. 20 

III. FOUR-SPIN EXCHANGES IN 
U(1)-SYMMETRIC SPIN-| LADDERS 

As additional four-spin exchanges to the U(l)- 
symmetric spin ladder (1), let us consider the following 
terms: 

Hvl = Mi ^(Sij ■ S 1>j+1 )(S 2 ,j ■ S 2 ,j+i), (6a) 

3 

Ji„ = Vrr y^(Sl,j ' S 2 ,j){SlJ + l ■ S 2 ,j+l), (6b) 

3 

H x = V x ^(Su • S 2 j +1 )(S 2tj ■ S w ). (6c) 

3 

The well known four-spin cyclic term satisfies V\i = V IT = 
— V x - 3 The bosonized spin operators (2) and symmetry 
arguments provide the following bosonized expressions of 
four-spin exchanges: 

/dx 
— Cncos(>/8^-) + -" , (7a) 
ClQ 

/dr 
— M 2 CVicos(V2^6L) 
a 

+C r2 cos(2\/27r0_) + C rZ cos(\/87rV>_ ) 
+C r4 cos(2\/8^0_) + • • • , (7b) 

/dr 
Af 2 C x icos(\/2^_) 
a 

+C x2 cos(2V2^6»_) + C x3 cos(\/8^>_) 

+ C x 4COs(2v / 87T(/>-) H , (7c) 



where C/i, C rq and C xq are nonuniversal constants, and 
we have written down only important terms. We note 
that (i) C11, C rq , and C xq with q = 1,2, and 4 are 
all positive, (ii) C rq is the same order as C xq when 
q = 1,2, and 4, and (iii) CV3 and C X 3 depend strongly 
on the magnetization M. In Eq. (7a), we have used 
the bosonized dimerization operator (— l) 3 S n j -S n j+i ~ 
cos(v / 47T'/ , n + 2irMj) + ■ ■ ■ } The 9 + vertex operators and 
4> + ones without oscillating factors are absent in Eq. (7) 
because the four-spin exchanges (6) conserve all the sym- 
metries of the ladder (1). It is remarkable that 7i rr and 
TL X include a new 9- term, cos(2v / 2~7r0_), which can 
make the locking position of 6>_ shift. From the calcula- 
tion in Eq. (7), we find that a four-spin exchange term 
with the form 

{S\,% ■ S 2 j)(Si ! k ■ S 2 ,i) (8) 

always contains cos(2\/27r#_). 21 The scaling dimension 
of cos(2\/27r#-) is 2/K^{k 2/K), and thus it becomes 
relevant when K- > 1. Equation (4) and the property of 
K suggest that the condition K_ > 1 is realized when the 
system is sufficiently close to the saturated state (K — > 1) 
and the rung couping is AF ( J± > 0). If we introduce 
perturbations such as further-neighbor exchanges for the 
intrachain direction, the parameter regime with K- > 1 
could expand. In addition, four-spin exchanges would 
vary the value of K_ slightly. 

Let us consider a spin ladder with four-spin terms V TI 
and Vx under the condition K_ > 1. In its bosonized 
Hamiltonian, the TLL of the (<j)+,9+) sector is still sta- 
bilized by symmetries of the ladder. On the other hand, 
the {(f)- , 9^ ) sector is described as 

^[4>-,0-] = Gaussian part (9a) 
+v_ { dx cos(V2^9) + cos(2V2^l9_), 

f 51 oc B$J ± + M 2 (C r iV rr - CxiVx) , . 

\ 92 oc C r2 V„ + C X2 V X ' [yD > 

where a ~ ag is the short-distance cut off, g\ y2 are dimcn- 
sionless coupling constants, d\ = 2 — \j(2KJ) and d 2 = 
2 — 2/K_: 1/(2K_) and 2/K_ are the scaling dimensions 
of cos(v / 27r#-) and cos(2\/27r^_), respectively. Here we 
have omitted irrelevant 0_ vertex terms. Equation (9a) 
is a well-known double sine-Gordon (DSG) model. 22 It 
is widely believed that the DSG model exhibits an Ising- 
type quantum phase transition. In fact, one easily finds 
that when g 2 is positively increased enough compared 
with |<7i|, the potential form of the DSG model changes 
from a single-well to a double-well type. The Ising transi- 
tion takes place just between the single-well (disordered) 
phase and the double-well (ordered) one. For the or- 
dered phase, the locking position of deviates from 
and ±y/ir/2. It indeed means the emergence of a fi- 
nite vector chirality (Vj) and the spontaneous breaking 
of the Z 2 rung-parity symmetry. Note that even in the 
chiral-ordered phase, the U(l) symmetry around the S z 



4 



axis (8+ — > 6 + + const) is conserved. It is expected that 
even via the Ising transition, there is no singular behavior 
in the magnetization curve since the DSG model in the 
(</>_, 6L) sector does not directly couple to the field H. 
Equation (9b) indicates that the system favors the vector 
chiral long-range order, for example, under the following 
conditions: 



ordered phase 
(vector chiral 
order in TLL) 



(1) V„ >J±>0 and V x ~ 0, 

(2) V x > J± > and V„ ~ 0, 

(3) V rr ~ V x > Jj_ > 0, 



(10) 




Ising 
fixed point 



disordered phase 
(field-induced TLL) 



Gaussian fixed point I £>i I 

FIG. 1: RG flow and phase diagram of the DSG model with 
K- > 1 around the Gaussian fixed point ( (71,32) = (0,0). 
Phrases in the parentheses stand for the phases of the cor- 
responding spin-| or spin-1 AF ladders with four-spin ex- 
changes and a Zeeman term: 31,2 are given by Eq. (9b) 
[Eq. (17)] inthespin-i [spin-1] case. Near the Gaussian point, 

the RG-invariant curve satisfies g 2 ~ \gi\ d2,dl . The RG anal- 
ysis of the DSG model has been discussed in Refs. 24 and 
25. 



etc. It also tells us that the cyclic exchange term with 
V rr = — V x possesses only a low possibility of causing 
the chiral order. From these arguments based on the 
bosonization, we conclude that the vector chiral long- 
range order emerges in the field-induced TLL, when a 
sufficiently strong four-spin exchange (8) is introduced 
in the ladder with K_ > 1. Furthermore, if K- > 1 
and M — (H < Hi) simultaneously hold in a U(l)- 
symmetric ladder, i.e., if K- > 1 holds when the ((j>+, 0+) 
sector has a gapped spectrum (although this situation is 
hard to occur within a realistic ladder system), it is also 
possible that a four-spin exchange (8) yields a gapped 
chiral-ordered state. 

In order to see the global phase structure of the DSG 
model around the trivial Gaussian fixed point (51,(72) = 
(0,0), the renormalization-group (RG) analysis is very 
useful. Using the operator product expansion tech- 
nique, 23 we obtain the following one-loop RG equation where C n are nonuniversal constants, and £i( 2 ) is the cor- 



for the DSG model: 24 ' 25 



dgi 
dL 

dm 

dL 



= 2- 



1 



2K 
2 



3l - 7T5152, 



2" J—) - ; <7<I- 



(11) 



relation length in the disordered (ordered) phase of the 
DSG model. The result (12) would be useful in numeri- 
cally detecting the chiral-ordered phase. Furthermore, if 
the chiral-ordered phase is realized in a ladder compound, 
these features of the vector chirality could be observed 
experimentally, in principle. 



where L is the scaling parameter: a — > ae L . 
The RG equations have two nontrivial fixed points 
(±\/ d\ g?2 /ti") d\ 1 7r) = (±31,32). These must correspond 
to the above-mentioned Ising transition fixed point al- 
though their values g\ 2 would not be reliable. Analyzing 
the RG equations around (51,32) = (0,0) and (zbgfjpl)) 
we can draw the RG flow and the ground-state phase 
diagram as in Fig. 1. As expected, it shows that if 92 
sufficiently grows, the vector chiral order appears. 

The vector chirality Vj must play the role of the order 
parameter near the Ising transition. Therefore applying 
known results of the two-dimensional Ising model 9 ' 26 we 
can predict several properties of the chirality Vj. For 
instance, the chirality increases as (Vj) ~ (32 — ffl) 1 ^ 8 
in the vicinity of the critical point (31,32)- l n addition, 
the chirality correlation function near the critical point 
is predicted to behave as follows: 



lim (VjV ) 



C ie -bV£i/ 
Co/bl 1 / 4 
C 3 (32 



3 2 c ) 1/4 + C^e-I^A 



(52 < 32) 

(32 = 3 2 c ) ,(12) 

171 (92 > g\) 



As we already mentioned in the Introduction, in spin- 
i ladder systems, there is a duality relation that has 
been studied in Ref. 3. Let us here apply the duality to 
our spin- i ladders with four-spin exchanges. The duality 
transformation is defined as Sij = ^(Tij +T2,j) +T] 



T 2j and S 2tj = !(Tij + T 2j ) - 7,., x 
tors T n j obey the same algebra as spin- 



T2J. New opera- 
| operators. This 
mapping therefore makes an arbitrary spin- i ladder with 
{S n j} change into a new dual ladder with {T nj }. It is 
remarkable that the duality transformation conserves the 
total spin in each rung (Sij + S2J = T\j + T2,j) and 
changes the vector chirality in each rung to a Neel-type 
moment (— 2S±jX S2J = Tij—Tsj)- On the other hand, 
recall that a spin- 1 ladder H* = Hi a d[</, + ^rr[Vr] + 
H x [V x ] is shown to possess the vector chiral long-range 
order near the saturation if the condition > 1 is satis- 
fied. Combining this prediction and the above properties 
of the duality, we can find that the dual model for 
has a rung Neel order (T* • — Tf j) 7^ in the vicinity of 

1). The Hamiltonian of the 



the saturation (Tf + T : 



2, j 
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dual ladder is represented as 

,J Vy 



Wdual = X! ' T « J + l + X! Tl J ' T2 -- 



+ (-. 



J Vx 



( 2 J + T") E( T M • T w)(T2,j ■ T W ) 



V, 



fields <j) n and # n : a U(l) rotation around the S z axis 
and the one-site translation respectively correspond to 
n -» #n + const and (</>„, 6>„) -> (<£„ - M^/n, 6 n + V^). 30 
These symmetry operations are very similar to those of 
the spin-i case. 

Substituting the formula (14) into the rung coupling 
of the spin-l ladder (1) and integrating out the massive- 
magnon part in its effective Hamiltonian, 32 we obtain a 
two-component boson theory that is the same type as 
Eq. (3). Namely, the (<f> + , 6 + ) sector again brings a TLL, 
and the (</>_ , ) sector takes a massive spectrum with 
pinning 6L to or ±^/tt/2. 33 The TLL state is strongly 
protected by symmetries of the ladder like the case of 



,j+i) the spin-i ladder. Formula (14) further tells us that the 



(13) 



Since this model has a strong rung-coupling term and 
strong four-spin ones, it cannot be analyzed within the 
weak-rung-coupling approach in this paper. 



IV. SPIN-1 LADDERS 

Let us now turn to the two- leg spin-l ladder (1) where 
spin-i operators are replaced with spin-l ones. One will 
encounter a scenario similar to the spin-^ case below. 

In the decoupled case with zero field, each AF chain 
has a Haldane gap, ~ 0.41J. When the magnetic field 
H exceeds it, a field- induced TLL appears due to the 
S z — 1 magnon condensation. For this TLL phase, the 
effective field theory has been established. 27-30 The ef- 
fective Hamiltonian consists of the Gaussian part, which 
corresponds to the TLL, plus the part of S z = and 
— 1 massive magnons. In contrast to the spin-i case, the 
TLL parameter K is larger than l. 31 The effective theory 
enables us to write the field theory form of spin operators 
as follows: 

S*j ~ M -a d x <j> n /^ + D l cos{s/4^<t> n -2-n Mj) 
+ {-l) j [D 2 <7 n cos(V^0n - nMj) 
+£) 3 (r ? „e^"+h.c.)] + --- , (14a) 



°n,j 



(-iy e- l ^ e ^ n [E 2 + E 3 cos(V4^„ - 2kMj)] 
+Eie -iV*e n cos (^ n _ nMj)(^ n + i£ fli „) 
+ •••, (14b) 



where D q and E q are nonuniversal constants, 



are the massless boson fields for the S z — 1 condensed 
magnon, the set (£„, n , c„, Mn) describes the S z = 
magnon sector, and rj n is the S z — — 1 magnon field. For 
the massive S z — magnon sector, (a n ) = and (/i n ) ^ 
hold. [For a more detailed explanation of Eq. (14), 
see Ref. 30.] If we consider only part of the TLL sector, 
the symmetries of the spin-l ladder can be expressed by 



vector chirality is written as 



V-j ~ ^1^2 sin(v27r#_) + 



(15) 



and leading parts of Af±j are written as a function of <j>±. 
These field theory results are indeed analogous to those 
of the spin-i case, Eq. (5). The formula (15) shows that 
(Vj) = holds in the standard spin-l ladder (1) with 

For the TLL phase of the spin-l ladder, let us add a 
four-spin interaction term. In this paper, we focus on a 
so-called biquadratic term 11 ' 12 



(16) 



that is one of realistic, familiar four-spin terms in spin- 
1 systems. Employing the field theory formula (14) 
and symmetry arguments, and then tracing out massive 
magnons, we obtain the bosonized form of Tib, which con- 
tains cos(2\/27r#-). As a result, the (0_,6*_) sector in- 
cluding effects of Hb could be described by a DSG model 
that has the same form as Eq. (9a). For the present case, 
the coupling constants are evaluated as 



.91 oc dJ±+M 2 C 2 V 
52 oc C 3 V + Ox (j) V 



O2 [ J ~j ] ■/_. 



(17) 



where C q are nonuniversal positive constants, and the 
second and third terms of g 2 originate from the trace- 
out procedure of massive magnons. As we already men- 
tioned, since the TLL parameter K of the spin-l AF 
chain is larger than 1, K- > 1 would hold in a wide 
parameter regime where cos(2v / 2~7T0-) is relevant. From 
Eq. (17) and Fig. 1, we can conclude that if g 2 > \gi\, 
namely, V > | J±\, under the condition K_ > 1, a vector 
chiral long-range order emerges like the spin-i case. In 
addition, Eq. (15) indicates that around the Ising tran- 
sition, the chirality correlation function in spin-l ladders 
behaves as Eq. (12). As one easily expects from Eqs. (8) 
and (14), besides the biquadratic term, other four-spin 
exchanges with the form (8) can also cause the vector 
chiral order. 
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V. PARITY-BREAKING PERTURBATIONS 

From Sec. II to Sec. IV, we have studied two-leg spin 
ladders with the Z 2 rung-parity symmetry. In this sec- 
tion, we shortly discuss effects of rung-parity-breaking 
perturbations. Let us focus on the following three kinds 
of realistic perturbations for spin- i ladders: 

Hudm = ^BufSijX^) 2 , (18a) 

3 

H S BM = ^(-l) J 'l> a (5ij x S 2J ) Z , (18b) 

j 

^chiral — ^ ^ -^pqr^p 

(S q x S r ), (18c) 

where p, q, and r in Eq. (18c) represent a site index on a 
spin ladder, {n,j}. The uniform DM term (18a) and the 
staggered one (18b) 13 often emerge in sufficiently low- 
crystal-symmetric magnets, and their origin is spin-orbit 
coupling. The third three-spin term (18c) 14,15 might be 
unfamiliar. It originates from virtual electron-hopping 
processes on a triangle plaquette {p, q, r} in a Mott- 
insulating ladder under an external magnetic field. The 
coupling constant A pqr oc sm(e$ pqr ), where e is the elec- 
tron charge and <fr pqr is the magnetic flux enclosed by 
the plaquette {p,q,r}. It therefore depends heavily on 
the direction and the strength of the magnetic field, and 
vanishes when the plaquette plane and the direction of 
the magnetic field is parallel or the magnetic field is ab- 
sent. (For a more detailed discussion on the three-spin 
term, see Rcfs. 14 and 15.) 

Formula (2) allows us to bosonize the DM terms as 
follows: 

/dx 
— £2 s in(V2^9_) + --- , (19a) 
a 

Hsbm « -D. ( — (-iyB%sm(V2^0-) + -- - .(19b) 
J «o 

This result is obvious from the bosonized expression of 
vector chirality, Eq. (5a). As we already discussed, the 
phase field 6>_ is locked at or ±y/w/2 in the spin ladder 
Wiad due to the potential cos(y/2n9-). Equation (19a) 
clearly shows that as soon as a uniform DM term is added 
to 7ii a d , the locking position of 9- immediately varies and 
the vector chirality takes a finite expectation value. In 
contrast, since a small staggered DM term is irrelevant 
due to the factor (— l)- 7 ', it hardly affects the low-energy 
physics of the ladder Hud- 

Similarly, one can bosonize the three-spin term (18c). 
It is readily found that Eq. (18c) also contains 
sin(\/27r#_), and is hence expected to cause a finite vec- 
tor chirality. However, one should note that besides the 
sine term, the three-spin term includes other relevant or 
marginal terms such as d x (f>±d x 6 T - Therefore effects of 
the three-spin term would not be simply argued. For in- 
stance, authors in Ref. 15 have obtained an unexpected 



prediction that for spin- \ easy-plane zigzag ladder, whose 
ground state has a vector chiral order, a three-spin term 
helps the chiral order vanish. The full bosonized form of 
Eq. (18c) and its effects on spin ladders depends on the 
coupling constant A pqr , namely, the underlying Mott- 
insulating ladder. In spin-1 or higher-spin ladder sys- 
tems, parity-breaking terms such as Eq. (18) would be 
also able to appear, and they induce a finite vector chi- 
rality. 

From the arguments in this section, we can expect that 
if a small parity-breaking term like Eq. (18a) is present 
in spin ladders with four-spin exchanges, it could be a 
trigger of the emergence of the four-spin-exchange- and 
magnetic-field-drivcn vector chiral order. Furthermore, 
such a term would help the chiral order survive even at 
a finite temperature. 



VI. CONCLUSIONS AND DISCUSSIONS 

We have studied effects of four-spin exchanges in two- 
leg spin- 1 and spin-1 AF ladders under a magnetic field, 
utilizing bosonization techniques. The magnetic field re- 
duces the spin SU(2) symmetry to the U(l) one, and as 
a result, the emergence of vector chirality (Vj) becomes 
possible. In Sees. Ill and IV, certain types of four-spin 
exchanges with the form (8) are indeed shown to yield the 
vector chiral long-range order (Vj) ^ with the Z 2 rung- 
parity symmetry spontaneously broken and the U(l) one 
unbroken, if the four-spin exchanges are large enough and 
the TLL parameter K- is larger than 1. Spin-1 ladders 
more easily satisfy the condition if_ > 1 rather than 
spin-i ones: for the spin-i case, K_ > 1 is predicted 
to hold only near the saturated state, while for the spin- 
1 case, it would hold in a wide region from the lower 
critical magnetic field to the upper one. An Ising-type 
phase transition takes place between the vector-chiral- 
ordered state and the standard spin liquid. Around the 
transition, the vector chirality correlation function fol- 
lows Eq. (12). The Ising transition would not accom- 
pany any singularity of the magnetization curve, because 
the (4>-, 0-) sector, which brings the transition, does not 
directly couple to the magnetic field H . It is notewor- 
thy that for spin-i ladders with a four-spin cyclic ex- 
change, the chiral order is hard to generate. For spin-i 
ladders with four-spin exchanges, we have also applied 
the duality transformation. 3 The dual spin ladder (13) 
is predicted to possess a Neel order for the rung direc- 
tion. In Sec. V, we have briefly considered effects of rung- 
parity-breaking perturbations on spin ladders. In partic- 
ular, we clearly find that a uniform DM interaction (18a) 
yields a finite expectation value of vector chirality in the 
ladder TL\ a< i even if it is small, while a staggered DM 
term (18b) is irrelevant in the same ladder. Moreover, 
it is expected that some types of these perturbations in- 
cluding Eq. (18a) could play the role of the trigger of the 
vector chiral order generated by four-spin exchanges and 
a magnetic field. 
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Our predictions in the spin- i and spin-1 cases suggest 
that for general spin- 5 AF ladders, a four-spin exchange 
such as Eq. (8) and a magnetic field can also induce a 
vector chiral order. Within the bosonization framework, 
it is difficult to quantitatively calculate critical values of 
the four-spin exchanges: V£, V£, and V c . Determining 
them accurately is an interesting future problem. 

Unfortunately, the predicted vector chiral long-range 
order has never been detected in experiments. However, 
since it is suggested 11 that effects of four-spin exchange 
cannot be negligible in some magnets and a magnetic 
field is one of the few things that we can control, there 
is a possibility that an experimental evidence of the chi- 
ral order is found. It is shown that a polarized neutron 
scattering experiment can detect a vector chiral order. 35 
In addition, if four-point spin correlation functions such 
as nonlinear susceptibilities are experimentally observed, 
they must contain some features of the vector chiral order 
or correlation [see Eq. (12)]. 

As is well known, theoretically, a typical origin of 
four-spin interactions is fourth- or higher-order electron- 
hopping processes in Mott insulators. 34 Hence the low- 
energy physics of Mott insulators with a large hopping 
matrix element would be described by a spin system with 
considerably large four-spin terms. Recently techniques 
in optical lattice have been greatly developed and they 
make it possible to construct several insulating states. 36 
Experiments in optical lattice therefore might succeed in 
realizing a vector-chiral-ordered state. It is also known 
that spin-phonon couplings can induce four-spin terms. 37 
Therefore magnets with such a strong coupling might be 
another candidate for a chiral-ordered state. 

It is shown in a previous work 38 that in spin-i lad- 
ders, a four-spin cyclic exchange can yield a M = 1/4 
plateau state with a Neel-type order when J± ^> J. In 



our weak-rung-coupling approach, such a plateau phase 
could also appear if a vertex operator cos(2v / 87T0+ + 
8wMj) = cos(2v / 87r</>+) in the (<j) + ,9 + ) sector becomes 
relevant, i.e., K + is smaller than 1/4. Supposing that 
K + < 1/4 and K- > 1 simultaneously hold in the 
weak-rung-coupling regime, we can expect that the vec- 
tor chiral order, caused by four-spin exchanges, still 
survives in the plateau state. In this plateau regime, 
there is the possibility that a subleading term of Vj, 
sin(v / 27T0_) cos(v / 87r^+ + AirMj), causes a staggered vec- 
tor chiral order. We, however, note that the conditions 
K + < 1/4 and K_ > 1 arc difficult to occur in real lad- 
ders. 

Spin- 5 ladders with four-spin exchanges possess an ex- 
actly solvable SU(4)-symmetric point, around which we 
can construct its effective field theory 4,8 The relationship 
between the SU(4)-point field theory and the bosonized 
one in the present paper has not been established well. 
Making it clear is an interesting open problem, and must 
contribute to a more sophisticated understanding of spin- 
\ ladders. 

We finally note that besides a magnetic field, other 
U(l)-symmetric terms such as XXZ anisotropy would 
help the emergence of the chiral order in AF spin ladders 
with four-spin exchanges. 
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